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Abstract We revisit the analogy suggested by Madelung between a non-relativistic 
time-dependent quantum particle, to a fluid system which is pseudo-barotropic, 
irrotational and inviscid. We first discuss the hydrodynamical properties of the 
Madelung description in general, and extract a pressure like term from the Bohm 
potential. We show that the existence of a pressure gradient force in the fluid 
description does not violate Ehrenfest’s theorem since its expectation value is 
zero. We also point out that incompressibility of the fluid implies conservation of 
density along a fluid parcel trajectory and in ID this immediately results in the 
non-spreading property of wave packets, as the sum of Bohm potential and an 
exterior potential must be either constant or linear in space. 

Next we relate to the hydrodynamic description a thermodynamic counterpart, 
taking the classical behavior of an adiabatic barotopric flow as a reference. We show 
that while the Bohm potential is not a positive definite quantity, as is expected 
from internal energy, its expectation value is proportional to the Fisher information 
whose integrand is positive definite. Moreover, this integrand is exactly equal to 
half of the square of the imaginary part of the momentum, as the integrand of the 
kinetic energy is equal to half of the square of the real part of the momentum. This 
suggests a relation between the Fisher information and the thermodynamic like 
internal energy of the Madelung fluid. Furthermore, it provides a physical linkage 
between the inverse of the Fisher information and the measure of disorder in 
quantum systems - in spontaneous adiabatic gas expansion the amount of disorder 
increases while the internal energy decreases. 
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1 Introduction 

A year after Erwin Schrodinger published his celebrated equation, Erwin Madelung 
showed (in 1927) that it can be written in a hydrodynamic form PQ. Madelung’s 
representation has a seemingly major disadvantage by transforming the single lin¬ 
ear Schrodinger equation into two nonlinear ones. Nonetheless, despite of its ad¬ 
ditional complexity, the hydrodynamic analogy provides important insights with 
regard to the Schrodinger equation [2j[3ll31[5l[6] . The Madelung equations (ME) 
describe a compressible fluid, and compressibility yields a linkage between hy¬ 
drodynamic and thermodynamic effects. The work done by the pressure gradient 
force to expand the flow transforms internal thermal microscopic kinetic energy to 
the macroscopic hydrodynamic kinetic energy of the flow. In this paper we wish 
to examine to what extent such a linkage can be made and what added value it 
provides in understanding quantum systems. 

ME can be obtained when taking the non-relativistic time dependent Schrodinger 
equation (TDSE) of a particle with mass m, under the presence of an external po¬ 
tential U(r, t): 




= 




V. 


(1) 


Assuming that the wave function 'I' is continuous and can be written in the polar 
form 

*(r,t) =Xr,f)e lS(r ’ t)/? \ (2) 

then together with the de Broglie guiding equation for the velocity 

u = VS (3) 


(where the tilde superscript represents hereafter a quantity per unit mass m, so 
that S = £-), the real part of the TDSE becomes the continuity equation 

In p = —V • u , (4) 

and the imaginary part becomes 

™ =H K + Q + U), (5) 

where K = u 2 /2 is the kinetic energy per unit mass and 


n 2 v 2 yp 


2 m 2 sfp 


( 6 ) 


is the Bohm potential per unit mass [7] . Generally (u • V)u = VA' + w x u, where 
uj = V x u is the vorticity, however for a potential flow in the form of (3), the flow 
is irrotational, i.e., lo = 0, and therefore when applying the nabla operator on (5) 
we obtain 


777U = -VQ(p) - VC/, 


( 7 ) 
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where + u • V, is the material (Lagrangian) time derivative of a fluid 

element along its trajectory. Equations (4) and (7) map the TDSE to a pseudo- 
barotropic, inviscid flow where (7) is its inviscid Navier-Stokes (i.e., Euler) equa¬ 
tion. As indicated by [8], in order for this formalism to be equivalent to the 
Schrodinger equation, one should also add a quantization condition. 

The paper is organized as follows. In Section 2 we examine the hydrodynamical 
properties of the Madelung equations. In Section 3 we relate them to their thermo¬ 
dynamical properties, and show how the Fisher information provides an analogy 
to the internal thermal energy of the flow. Discussion of the results appears in 
Section 4. 


2 Hydrodynamic properties of the Madelung equations 


2.1 The continuity equation 


It is somewhat surprising that we can describe the quantum state of a single parti¬ 
cle in terms of a fluid whose mass density pfi u id(r , t) = 5M/5V is the probability 
density p(r, t) of the wave function to find the particle m in location r at time t 
(5M is an infinitesimal “fluid mass” occupying an infinitesimal volume 5V, where 
M is non-dimensional since it represents probability). The continuity equation (4) 
is usually represented in its Eulerian form 


dp 

dt 


-v J, 


(8) 


where J = pu is the flow mass flux. In its Lagrangian representation however, 
(4) simply reflects the statement that a “fluid parcel” conserves its mass 5M as 
it moves with velocity u, i.e., -j=^( 5M ) = 0. (4) is then obtained when noting 
that V • u = -jj-j: ln(<5y) is the flow compressibility term. Hence, when the flow is 
incompressible -^-p = 0, and that is to say that the density is conserved along 
a “fluid parcel” trajectory. This is indeed the case when the wave function is 
represented by a single plane wave since then u = S,k, where k is the wavenumber 
vector. Nonetheless, interference between two plane waves or more yields V • u / 0 
in general, and thus the probability to find a particle along a trajectory constructed 
from (3) varies along the trajectory itself. In ID, incompressibility implies that 
u is not a function of space, and therefore p is non-spreading whether or not the 
flow is accelerating. 


2.2 Relating the Bohm potential to a pseudo barotropic pressure 

For an inviscid flow, in the presence of an exterior potential U, the Euler equation 
(Newton’s second law essentially) reads: 

TTt U= -)vp-w, (9) 

where the pressure P is a known thermodynamic property. Furthermore, if the 
flow is barotropic, that is the pressure P is a function of density only, the pressure 



4 


Eyal Heifetz, Eliahu Cohen 


gradient force (PGF) can be written as a perfect gradient 

- -VP(p) = -VQ(p), (10) 

P 

where Q = f is defined up to an unspecified time dependent gauge function. 
Hence, taking then the curl of (9) we obtain 

| = Vx(ux,), (11) 

therefore, if ui(t = 0) = 0, the vorticity remains zero at all times and thus the 
velocity field can be represented as a potential flow of the form of (3). In regions 
where p = 0, the vorticity can be non-zero in principle, as in superfluid singulari¬ 
ties, however we refrain here from discussing such cases. 

It is worth noting that Euler had derived equation (9) in 1757 more than 
a century before the thermodynamic kinetic theory of gases was established by 
Maxwell and Boltzamnn in 1871. Hence, the knowledge that the pressure results 
from the aggregated effect of microscopic random motions of molecules or atoms 
impacting each other, was not necessary to Euler to formulate it as a macroscopic 
force per unit area that is somehow associated with the intrinsic properties of the 
fluid. In a sense, this is what we try to do next with respect to the Madelung 
equation (7). Here, the starting point is the known function of the Bohm potential 
Q(p), hence we can solve (10) for P to obtain: 

p = n(p) + m (12) 

where 

< i3) 

and f(t) is a time dependent gauge function. Here we used the identity 

—- = V 1 2 lna+ (Vina) 2 , (14) 

a 

to extract 77 from Q in (10). As 77 depends on spatial derivatives of the density, 
it cannot be defined as a proper thermodynamic pressure, hence a possible ap¬ 
propriate name for it could be a pseudo pressure or a pressure like term Q. An 
example that makes sense of (13) is a normal probability density function, as it 
is straightforward to see that 77 oc p, which is true for an ideal isothermal perfect 
gas. 

In the absence of an external potential U, it is surprising that the quantum 
state of a free particle is mapped into a fluid that can be accelerated by a force, 

1 Indeed, more complex interpretations have been suggested with regards to the quantum 

pressure 131151 . however the proposed pressure like term still seems to us the most straight¬ 
forward one as we are seeking for the simplest analogy. A second rank tensor of the Cauchy 
stress tensor, embedding shear terms, as well as the introduction of turbulence to the ME, or 
the representation of it in terms of nonlinear diffusion are interesting interpretations, however 
they involve an amount of complexity that we try to avoid when considering the dynamics of 
a single quantum particle. 
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even if this force is the PGF which is internal to the flow. However, defining the 
acceleration as a = its expectation value can be found to be 


(a) = J ^*adPdV = J padV = — J 


PdA 


(15) 


(where the integrations are taken over the whole domain V). Thus, if p is bounded 
within the domain and vanishes at the domain boundary A , then the net expec¬ 
tation value of the acceleration is zero and Ehrenfest’s theorem is not violated. 


2.3 Flow incompressibility and non-spreading wave packets 

As pointed out in Section (2.1), when the flow is incompressible (V • u = 0), 
density is conserved along “fluid parcel” trajectories. It is worth noting that this 
statement is equivalent to the one that a wave packet following such trajectory is 
non-spreading. Furthermore, since ^ x u + V(u 2 /2), and the flow is 

irrotational, (7) becomes 

^ = -v[k + Q + u). (16) 

Specifically, for ID flow, incompressibility (§^ = 0) implies that u can be only a 
function of time and (16) becomes 

S = < 17 > 

posing the constraint on the potentials: Q + U = a(t)x + b(t). Among the potentials 
mentioned in the literature which satisfy this non spreading condition are: i ) the 
free Airy wave packet (yT> oc Ai(x) & U = 0) [9]; ii) the gravitational “quantum 
bouncer” (y/p oc Ai{x) & U oc x) TO]; in) the ground state of the harmonic 
oscillator (y/p oc e~ x & U oc x 2 ) [T0| . 


3 Thermodynamic like properties of the Madelung equations 

In this section we examine to what extent one may associate thermodynamic like 
properties with the Madelung description, in reference to the thermodynamic of 
classical barotorpic conservative (adiabatic) flows. 


3.1 Energy conservation 

Consider a classical, adiabatic (entropy conserving, hence inviscid), barotopric 
flow. In the absence of heat exchange, the first law of thermodynamics reads, 
dl = —PdV, where I is the thermal internal energy, representing the macroscopic 
aggregated effect of the microscopic random thermal fluctuations. Hence, during 
an adiabatic process, compression of a fluid parcel by its surroundings performs 
work that increases its internal energy. If, however, the flow is incompressible, the 
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internal energy remains unchanged. When following materially a fluid parcel in 
motion, the adiabatic first law is transformed into 

p R-I = -pR-\n{5V) = -PX/-u. (18) 

Multiplying (9) by pu and combining it with (18), then for a time independent 
external potential U, we obtain 

p-^- t (K+~I + u) =-V-(uP). (19) 

Thus, the total energy of a fluid parcel per unit mass 

Eci = {K + I + U), (20) 

is not materially conserved simply because the fluid parcel is not an isolated system 
(the subscript Cl denotes classical fluid). The surrounding pressure can change the 
parcel internal energy by compression and the pressure gradient can accelerate the 
parcel and hence change its kinetic energy. Nonetheless, the overall total energy 
of the fluid is conserved in the domain averaged sense. We can use (4) to obtain 
that for any scalar field a, p^ppp = J ^(pa ) + V • (pua), and therefore if all fluxes 
vanish on the domain boundaries we obtain 

w (M = ml piiadV =°- (21) 

On the other hand, direct calculation of the energy expectation value from the 
Schrodinger equation (1), yields the conservation of the total energy (assuming as 
well that all fluxes vanish on the domain boundaries): 

(*«■) - k / ' r ^ mv =<*+«+ 6 ) • (22) 


where Q is the Bohm potential (per unit mass) given by (6), and the subscript Qu 
refers to a quantum like fluid. Comparing (20) with (22) suggests therefore that 
(Q'j — A similar argument, arising from a different perspective, has been 

recently suggested in m The Bohm potential by itself, however, is not a positive- 
definite term, as is expected from internal energy. Nonetheless its expectation value 
yields, after integration by parts, 




Av(lnp) 



1 

2 



(23) 


where FI = f i(V p) 2 dV = f ^^(lnp)] 2 ^^ is the Fisher information, in the form 
presented in dim For instance for normal distribution the Fisher information 
is inversely proportion to its variance, and according to the Cramer-Rao bound 
the general variance of any unbiased estimator (like the location of a particle r at 
time t) is always larger (or equal for normal distribution) to 1 /FI [13] . Therefore 
we can choose to define the positive definite internal energy per unit mass as 








Title Suppressed Due to Excessive Length 


7 


Furthermore, this choice suggests a relation between the quantum thermal energy 
and the imaginary part of the momentum field. Writing the momentum in the 
form of 

pP = [—ihV lntf']'?' = p P, (25) 

we can define the complex velocity field as 

v = p = v r + ivi (26) 

where 

v r = u = VS, Vi = ——V(lnp). (27) 

2m 

This partition has been suggested, by mm in different contexts. Imaginary mo¬ 
mentum is a somewhat strange concept, however here we simply note that in the 
same sense that under the de Broglie guiding equation (3), the macroscopic ki¬ 
netic energy of the flow is K = ^v 2 , under (24) / = |v 2 , suggesting that within 
this analogy, the microscopic thermal motion is represented by the imaginary mo¬ 
mentum. Moreover, the expectation value of random fluctuations must be zero by 
definition, and indeed the expectation value of each component of is propor¬ 
tional to the Fisher score which is zero. 

As pointed out by mm the inverse of the Fisher information can be used 
to measure the degree of disorder in the system and it tends to increase with 
time. In this sense, the Fisher information measures the “narrowness” of the po¬ 
sition distribution around the actual position of the particle m ■ The relation of 
the internal energy of the Madelung flow to the Fisher information provides a 
physical interpretation of this statement. Consider a gas expanding spontaneously 
in an adiabatic process. The degree of disorder increases and according to the 
first law of thermodynamics (18) the internal energy decreases. The fact that the 
Fisher information may provide better estimation to disorder than the negative 
entropy (negentropy), may result from the fact that during adiabatic expansion 
the thermodynamic entropy remains unchanged. Furthermore, in the absence of 
an external potential, (^K + J^> is conserved and if <^/^> decreases during expansion 
the kinetic energy and momentum become more pronounced. This may suggest a 
different angle to the relation between the Fisher information and the Heisenberg 
uncertainty principles )16| . 


3.2 The Bohm potential and pseudo enthalpy 

Equation (10) can be rewritten in the form 


VQ = -VP= -5 Vp + V - 


(28) 


and for classical barotropic adiabatic flows we can substitute the continuity equa¬ 
tion (4) with the first law of thermodynamics (18) to obtain 


D_ ~ _ P_D_ 
Dt p 2 Dt^ 


±10.0.)-Sfi. 


(29) 
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Hence for classical barotropic adaibatic flows, (29) implies that 


Q — (I ~b —) "b Bc{ti — Ent BcCti , (30) 

L V p) id 

where Ent = is the enthalpy per unit mass and Be(t) is a time dependent 

gauge function that is equal to the Bernoulli potential (as will be seen in the next 
subsection). On the other hand, for the quantum flow, (6), (13) and (24) indicate 
that the Bohm potential satisfies 


Q 



(31) 


This modified form of quantum enthalpy does not violate the equality ^Qy = y/y, 

since it is straightforward to show that f IIdV = 2 Hence, although II{p) is 
not a positive definite quantity, its domain averaged value is always positive. 


3.3 The Bernoulli and Hamilton-Jacobi equations 


For completeness we discuss the relation between the Bernoulli and Hamilton- 
Jacobi equations in classical and quantum systems. If the flow is potential of the 
form of (3) we can obtain from (7) the time dependent Bernoulli equation 


V-Be(f) = 0, Be(t) 
so that together with (30) 


— + (K + Q + U), 


dS 

dt 


(K + Ent + t/) 



(32) 


(33) 


where Hci differs from Eci in (20) by the term of -y. Moreover, since the thermo¬ 
dynamic pressure is positive definite it is clear that ^77^} ^ (^)] ' With this 

alert in mind, (34) can be regarded as the fluid mechanics version of the Hamilton- 
Jacobi equation. Since the velocity potential is only a function of time and space, 
S = S( r, t) (and not of the momentum p), the material derivative = yy +2 K, 

is the total time derivative in the phase space of (r, p). Substituting in (33) we 
obtain 

^ = K - (Ent + t/) =L => S = j LDt, (34) 

where L and the velocity potential S may serve respectively, as the Lagrangian 
and action of this system (where the material integration in time follows the fluid 
parcel). Requiring the action to be stationary by Hamilton’s principle, we obtain 
that (7) becomes the Euler-Lagrange set of equations. To the best of our knowledge 
this relation between the Bernoulli and the Hamilton-Jacobi equations has not 
been acknowledged in the literature. 

In the quantum realm the starting point is equation (5), thus Be{t )=0 and 
Hqu = (K — I + -y- + U). But unlike the classical counterpart, here indeed 
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4 Discussion 

The Madelung equations suggests an alternative appealing description to the 
Schrodinger equation by associating the behavior of a non relativistic quantum 
particle to the dynamic of a fluid. This provides an additional important intuition, 
although that in principle the nonlinear fluid dynamic equations are much more 
complex than the single linear Schrodinger one. Any compressible flow exchanges 
kinetic energy with thermal internal energy as both result from the motion of 
the particles composing the flow. The Madelung flow is compressible, however the 
Madelung equations do not provide an explicit description of the involved thermo¬ 
dynamical processes. Moreover, the Madelung flow is a (pseudo) barotropic, poten¬ 
tial (and hence irrotational), inviscid flow. The thermodynamics of classical flows 
which encompasses such properties, is simple, elegant and intuitive. Therefore, 
our motivation here was to examine to what extent we can relate to the Madelung 
flow, thermodynamic like properties. Throughout the text we repeat using terms 
such as “like” and “pseudo”, in order to keep in mind that the Madelung fluid is 
a pure analogy. As pointed out in section 2.1 the “fluid density” is the probability 
density function to find a single quantum particle in location r at time t. 

In the first part of this paper we revisited the Madelung hydrodynamic analogy 
to the Schrodinger description in the presence of an exterior potential. The aim 
was to present this alternative, and somewhat surprising view, in a way that is 
appealing both to quantum physicists and classical fluid dynamicists. Furthermore, 
a few examples have been briefly discussed (with respect to ID non spreading wave 
packets) in order to exemplify how the hydrodynamic perspective can be intuitive 
and helpful to understand the behavior of quantum systems. 

In the second part we suggested a complementary, thermodynamic descrip¬ 
tion to the Madelung equations, where the classical thermodynamics is taken as a 
reference. We find that conservation of energy, implied from the Schrodinger equa¬ 
tion, suggests that the expectation value of the Bohm potential is proportional 
to the expectation value of the internal energy in classical fluids. However, the 
Bohm potential is not a a positive definite function of density, as is expected from 
the internal energy of barotropic fluids. Nevertheless, since the expectation value 
of the Bohm potential is proportional to the Fisher information, and the latter’s 
integrand is a positive definite function of density, we related it to an equivalent 
thermal internal energy. This provides a direct physical link between the decrease 
in Fisher information and the amount of disorder increase spontaneously with 
time. According to the Cramer-Rao bound the variance of the Fisher informa¬ 
tion is at least inversely proportional to the variance of the probability density 
function of the particle location, and as disorder increases with time the Fisher 
information decreases and the probability variance increases. In the compressible 
adiabatic fluid analogy, spontaneous adiabatic expansion of fluids decreases their 
internal energy and increases the amount of disorder in the fluids. Since in adi¬ 
abatic processes the thermodynamic entropy remains constant by definition, this 
analogy may strengthen the claims that the (inverse of the) Fisher information is 
a more suitable measure of disorder in quantum systems than entropy. 

Moreover, the quantum momentum operator is generally a complex function, 
however the de Broglie guiding equation makes use only of its real component. 
Nonetheless, if the amplitude of the density is space dependent the momentum 
has an imaginary counterpart. We find it elegant that as the macroscopic kinetic 




10 


Eyal Heifetz, Eliahu Cohen 


energy of the quantum fluid is proportional to the square of the real part of the 
momentum, the internal energy is proportional to the square of the imaginary 
one. Furthermore, the expectation value of components of the imaginary part of 
the momentum is zero (proportional to the Fisher score), which is expected for 
random thermal fluctuations. This suggests that within the fluid analogy the real 
part of the momentum represents the domain averaged momentum of the fluid 
and the fluctuation from this mean is represented by the imaginary part. 

The classical thermo-hydrodynamics of a barotropic conservative fluid cannot 
be mapped however, as is, to the Madelung fluid. This results from the different 
starting points of the two systems. In the classical form the pressure is a well 
defined thermodynamic property that is a function of density itself (and not of 
spatial gradients of density). Then the momentum dynamics is described by the 
Euler equation where the pressure gradient force can be introduced as a perfect 
gradient of a function that is density dependent but can be defined up to some ad¬ 
ditional time dependent gauge function. Furthermore, an internal thermal energy 
exists and obeys the first law of thermodynamics and this implies that the density 
dependent function of the perfect gradient is the enthalpy of the system. Domain 
energy conservation is obtained when transforming the momentum equation into 
mechanical energy one and augmenting it with the first law of thermodynam¬ 
ics and the continuity equation. One can also obtain the Bernoulli (which is the 
Hamilton-Jacobi) equation from the Euler equation up to the undefined Bernoulli 
potential time dependent function. 

In the quantum case, the real part of the Schrodinger equation is mapped into 
the classical continuity equation, however its imaginary part is mapped into a 
Bernoulli like equation with zero gauge function and a well defined density de¬ 
pendent Bohrn potential which is not the Enthalpy of the system. Then the Euler 
equation for barotorpic flow is obtained when the nabla operator is applied on the 
Bernoulli like equation. Hence, although the Euler equation in the two systems 
looks the same, it is obtained in two fundamental different ways. We can then 
extract a pressure like form from the Bohm potential but this pseudo barotropic 
pressure has a function that depends on the density gradient rather then the den¬ 
sity itself, and moreover can be defined only up to a time dependent gauge function. 
Furthermore, no internal energy is defined a-priori from thermodynamic consider¬ 
ations, and the first law of thermodynamics is not a postulate of the system. On 
the other hand, the domain averaged conservation of the total energy is obtained 
directly from the Hamiltonian of the Schrodinger equation. Although it looks like 
the Bohm potential plays the role of the internal energy in the Hamiltonian it is 
not a positive definite quantity, as opposed to the Fisher integrand that provides 
the same expectation value as the Bohm potential. 

Overall, the suggested thermodynamic description of the Madelung flow is 
formed here using a top-down approach in the sense that we do not build the 
macroscopic thermodynamic parameters bottom-up from a comprehensive statis¬ 
tical microscopic theory. This is somewhat similar to the approach taken by Euler 
when deriving the simplified momentum equation of fluid flows before the ther¬ 
modynamic kinetic theory of gases has been developed. 

In future works it could be interesting to relate this suggested formalism to 
Zitterbewegung in a similar manner to m and also to the zero point field in QED 
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